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' Abstract 

^ ^ We consider the dynamics of a charged particle interacting with background electromagnetic 

field under the influence of linearized gravitational waves in the long wave-length and low- 
^ I velocity limit. Following the prescription in ,5j, the system is quantized and the Hamiltonian 

is then solved by using standard algebraic iterative methods. The solution is in conformity 
with the classical analysis and shows the possibility of tuning the frequency by changing the 
magnetic field to set up resonance. 

\ With the development of various ground based gravitational wave (GW) detectors like LIGO 
ff^ ■ [1], the possibility of direct detection of GW(s) with a strain sensitivity of the order of /i ~ 

, 10~^^/V Hz or better in the frequency range between 100-1000 Hz is expected in near future 

[2]. With such a small response of matter to a passing GW, more often than not, a quantum 
mechanical treatment of this interaction is desirable, if not essential [3ll3|. Quantum dynamics 
. of a test particle and a harmonic oscillator under the influence of linearized GW(s) has been 
investigated earlier [5] and recently we have generalized it in a noncommutative (NC) [6] setting 
in a couple of papers [71 18| where we treat the theory to first order in the NC parameter. These 

■ studies have revealed the possibility of identifying the NC nature of spacetime (or at least put 
^ , stringent upper-bounds on the NC parameters [Ql [10} W\\ [T2| [T3| I14j ) through GW detection 

data. 

Investigations related to the coupling of a background electromagnetic wave with a GW 
and more specifically, their conjugate effect on test matter, however, has been missing in the 
literature and is of foremost interest in its own right. Furthermore, in recent years there is a 
resurgence of interest in this age-old problem of the coupling of electromagnetic field to GW since 
experimental investigations are under way with ever increasing precision both in the laboratory 
or by astronomical observation [151 116j . Also, analytic treatement and numerical simulation of 
the interaction of a GW with a strongly magnetized plasma have shown that for strong magnetic 
fields (~ 10^^ Gauss), the GW excites electromagnetic plasma waves where the energy absorbed 
from the GW by the electromagnetic oscillations is comparable to the energies emitted in the 
most energetic astrophysical events, such as giant fiares on magnetars and possibly even short 
gamma ray bursts |17j . 

In this paper, we therefore construct the quantum mechanics of a charged test particle 
gyrating in the presence of a constant background magnetic field along the z-axis, while lin- 
early polarized GW(s) parallel to the magnetic field is passing. This system has been analyzed 
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classically by several authors [T^ where one solves the geodesic equation in a linearized GW 
background. The problem is worth studying because of its inherent importance in various as- 
trophysical contexts mentioned above. A quantum mechanical formulation is essential since this 
would allow us to compare our results with the classical results and give us deeper insight. Also, 
since GW couples to electromagnetic waves in plasma mainly by generating electric currents 
inside the plasma by perturbing the charged particle trajectories [19], the quantum mechan- 
ics of the charged particle will also be relevant in that context. Last but not the least, the 
quantum mechanical formulation would enable us to apply the established framework of NC 
electrodynamics [20l EH [22l [23] and gravity [Ml [SH [26] to elevate the present study at the NC 
level. 

In our earlier work [Tj |8], we have presented the argument |5j of working back from the 
geodesic deviation equation in the proper detector frame to the corresponding Hamiltonian. 
Following the same path, the Hamiltonian in the present case becomes 



H 



1 

2m 



Pj 



-Aj + mVl^Xk 



a) 



where the notations have their usual meaning. If the GW is propagating along the z-axis, due to 
its transverse nature, T^ok will have non-zero components only in the x — y plane. We therefore, 
confine our attention to that plane. Assuming the symmetric gauge Aj — ^ 
the traceless property of the GW [5], the Hamiltonian ([T|) takes the form 



-nejixi and using 
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2m 



2 2 

-mzu Xi 



wL + -hjkXkPj 



eB 
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eijhjkXiXk . 
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The first two terms represent an ordinary harmonic oscillator with the cyclotron frequency 
■07 = the third is the Zeeman term with L = eijXipj, the fourth term which is linear in 
the afhne connections, reveals the effect of the passing GW and finally the last term shows the 
coupling between the GW and the background magnetic fielcill- Note that, in the symmetric 
gauge we have obtained a Hamiltonian similar to that of a harmonic oscillator along with 
other terms which can be treated at different orders of perturbation. According to the classical 
treatment when the GW propagates parallel to the magnetic field with a frequency w' , the 
coupling of a gyrating particle with the GW becomes very strong at twice the cyclotron frequency 
zu' = 1w |18j . We shall try to find whether the quantum mechanical model also show such 
resonant behaviour. 

Defining raising and lowering operators in terms of the oscillator frequency w 



Pj 



2mw 



1/2 



(aj + a]) 



/ Timw 



1/2 



(3) 
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we write the Hamiltonian 



as 



H 



+ 1 ) + ihweija\aj 



eijhjk [aiak - aja|. 



+ 



'-rhjk ( 



(5) 



Since we are dealing with linearized gravity, a term quadratic in F has been neglected in eq.Q. 
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Working in the Heisenberg representation, the time evolution of aj(t) is given by 

dai{t) . 1. I i . / I 

— ^— = -iwai + weijUj + -riija'j + -eijhjk [a'f, + Ofc 

+\^kjhji {a\ + ak) (6) 

and that of a\{t) is the complex conjugate (c.c) of the above equation. Next, noting that the 
raising and lowering operators must satisfy the commutation relations 



aj{t),al{t) 



Ojk 



[aj{t),ak{t)] = 0=[a]{t),al{t)\ (7) 

we write them in terms of aj{0) and aj(0), the free operators at time t = 0, by the time-dependent 
Bogoliubov transformations 

aj{t) = Ujk{t)ak{0) +Vjk{t)al{0) 

a]{t) = al{0)ukj{t) + ak{0)vkj{t) (8) 

where the bar denotes the c.c and Ujk and Vjk are the generalised Bogoliubov coefficients. They 
are 2x2 complex matrices which, due to eq. ([7|), must satisfy uv'^ = u^v , uv) —vv'^ = /, written 
in matrix form where T denotes transpose, f denotes c.c transpose and / is the identity matrix. 
Since aj{t = 0) = aj(0), Ujk{t) and Vjk{t) have the boundary conditions 

Ujk{Q) = I , t;,fc(0)=0. (9) 

Then, from eq.® and its c.c, we get the following equations of motions in terms oi C, = u — v'^ 
and ^ = u + v'^: 



= -iwQi + weijE^ji + - 



iwC,il + roey^j7 + -hijiji (10) 



^^^^ ■ t , r- '^U r 

— - = -lWt,il + ^(-ijQjl - ^'T'ijQjl 
1 . 2, ■ 

+ -^ijhjkikl — -jhij€jkikl ■ (11) 

We shall solve the above eq(s) (fTO | [TTI) for the special case of linearly polarized GW(s). 

In the two-dimensional plane, the GW, which is a 2 x 2 matrix hjk, is most conveniently written 

in terms of the Pauli spin matrices as 

hjk (t) = 2/(t) (excrk + e+cj-fc) ; = ei, £+ = es (12) 



where 2/(t) is the amplitude of the GW and Ex (t) and e+ (t) represent the two possible po- 
larization states of the GW and satisfy the condition + = 1 for all t. In case of linearly 
polarized GW(s) however, the polarization states are independent of time and f{t) is arbitrary. 
To set a suitable boundary condition we shall assume that the GW hits the particle at t = so 
that 

f{t) = , for t < 0. (13) 

We now move on to solve eq(s) (llOl llip by noting that any 2x2 complex matrix can be written 
as a linear combination of the Pauli spin matrices and identity matrix. Hence we make the 
ansatz : 

Cjk{t) = AIjk + Bia]j^ + B2a% + B^a% (14) 
i,k {t) = Cljk + Dia]k + D2a% + Dsa% . (15) 
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Substituting for hjf^, (j^ and from eg (s) [HI in eg (s) PU] [TT ]) and comparing the coeffi- 
cients of / and fj-matrices, we get a set of first order differential eguations for A, Bi, B2, -B3, C, Di, D2,D^: 



A = 


-ivjC + iwB2 - f {eiBi + £3^3) - 2if {e^Di 


- eiDs) 


Bi = 


-iwDi -wB^-j {eiA - ^£3^2) + 2/ {eiD2 


- iesC) 


B2 = 


-iwD2 + iwA - if (£3^1 - £1^3) - 2/ {siDi 


+ £303) 


B3 = 


-iwD^ + TuBi-f {esA + ieiB2) + 2/ (£3£'2 


+ ieiC) 


C = 


-iwA + iwD2 + / {eiDi + £31)3) 




Di = 


—iwBi — wD^ + / (£iC — 18302) 




D2 = 


—iwB2 + iwC + if (£3-Di — £il?3) 




D3 = 


-i'wB3 + wDi + f{e3C + ieiD2) ■ 





(16) 

Noting tliat \f{t)\ « 1, tlie above set of eguations can be solved iteratively about its f{t) = 
solution by applying the appropriate boundary conditions ([H [T3|) . We therefore obtain to first 
order in the GW amplitude 



A{t) 
Bi{t) 



B3it) 

B2{t) 
Di{t) 



C{t) = e"*^* cos(rot) 

-/(t)e"*'^*{(£i + 2i£3) cos tut + (£3 - 2i£i) sin rot} 

+3ro (£3 - iei) f dt'f{t')e-^'^^' 
Jo 

-/(t)e~*^* {(£3 - 2z£i) cos rot - (£1 + 2i£3) sin rot} 

+3ro (£1 + i£3) /* dt'f{t')e-^'^^' 
Jo 

-ie-*"^* sin(rot) 



f{t)e (£1 cosrot + £3 sinrot) - ro (£3 - i£i) / dt f{t )e 



J./ f t -2ivjt' 



D3{t) = /(t)e~*^*(£3Cosrot-£isinrot)+ro(£i + i£3) / dt' f{t')e~^'^^' 

Jo 



(17) 
(18) 

(19) 
(20) 
(21) 

(22) 



Computing the above expressions for A, Bi, B2, B^, C, Di, D2, assuming a monochromatic 
GW waveform /(t) = /oe*^ * oscillating with frequency w' finally lead (by carrying out the 
procedure in [7]^]) to the following general expression of {Xi{t)) 



{Xi (t)) 



2h 



-Re[ai (t)] = Xi{0) cos^ rot + P,{of-^^^ + ^2(0) 
mro 2?7iro 



sin 2rot 



+ i^2(0) 



sin rot 



mw 
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mro 



fo 



F + £3 cos ro't + G 



(23) 



where {Rq) = (Xi(0), ^2(0)) and (Pq) = (-Pi(O), -P2(0)) are the initial position and momentum 
expectation values when the GW just hits the system at time t = and 



F(t) 
G(t) 



ro , ro , 

■ sin (2ro + ro ) t + 7: j sin (2ro — ro ) t 



2ro + ro' 
2ro2 



4ro2 - ro'2 4(2ro + w 
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- cos (2ro - ro') t . (24) 



4 



The above result is very similar to the one obtained earlier in [5] in noncommutative two- 
dimensional space and clearly shows that some of the oscillatory terms representing the response 
of the particle in the background electromagnetic field to the GW, will exhibit resonance when 
the condition w' = 2w is satisfied. This is indeed reassuring that the condition of resonance 
obtained through our quantum mechanical analysis is in conformity with the classical analysis 
|18j . Hence, one can in principle detect the GW and find its frequency w' by tuning the harmonic 
oscillator frequency w (which in turn depends on the applied magnetic field B) till resonance is 
achieved. Realistic scenarios can however be obtained if various forms of periodic GW signals 
with more than one frequency are used to carry out similar computations. We hope to report 
this in future. 
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